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Abstract 

In a diquark-diquark-antiquark picture of the pentaquark we study the decay — > K + n within 
the framework of QCD sum rules. After evaluation of the relevant three-point function, we extract 
the coupling gQ n K which is directly related to the pentaquark width. Restricting the decay dia- 
grams to those with color exchange between the meson-like and baryon-like clusters reduces the 
coupling constant by a factor of four. Whereas a small decay width might be possible for a positive 
parity pentaquark, it seems difficult to explain the measured width for a pentaquark with negative 
parity. 
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I. INTRODUCTION 



The experimental evidence for pentaquarks is now considered to be weak. However the 
strong statement that they do not exist can only be made after tests, which are presently 
being carried on [j| . Even if this state turns out to be unreal, its study generated many new 
theoretical ideas which can be valuable in the future Q|. 

One of the most puzzling characteristics of the pentaquark is its extremely small width 
(much) below 10 MeV which poses a serious challenge to all theoretical models. Many 
explanations for this narrow width have been advanced Q. In this work we calculate the 
pentaquark decay width within the framework of QCD sum rules (QCDSR) Several 
sum rule calculations have been performed for the mass of the pentaquark containing a 
strange quark These calculations are based on two-point functions with 

different interpolating currents. Surprisingly, all these determinations give similar masses 
with reasonable values. A common problem of all determinations is the large continuum 
contribution which has its origin in the high dimension of the interpolating currents and 
results in a large dependence on the continuum threshold. Another problem is the irregular 
behavior of the operator product expansion (OPE), which is dominated by higher dimension 
operators and not by the perturbative term as it should be. 

Here we present a sum rule determination for the decay width based on a three-point 
function for the decay O — > nK + . In this way we can extract the coupling ge n K which is 
directly related to the pentaquark width. To describe the pentaquark we use the diquark- 
diquark-antiquark model with one scalar and one pseudoscalar diquark in a relative S-wave. 

In (2I Q] it has been argued that such a small decay width can only be explained if 
the pentaquark is a genuine 5-quark state, i.e., it contains no color singlet meson-baryon 
contributions and thus color exchange is necessary for the decay. The analysis presented 
both in P| and in is only qualitative. The narrowness of the pentaquark width can 
then be attributed to the non-trivial color structure of the pentaquark which requires the 
exchange of, at least, one gluon. In this work we will also test quantitatively the hypothesis 
put forward in (t| ^| to see whether this mechanism is sufficient to explain the small width. 
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II. CORRELATION FUNCTION 



The investigation of the pentaquark decay width requires a three-point function which 
we define as 

r(p,p«) = Jd 4 xd 4 ye- l(iy e ip,x r{x,y), 

T(x,y) = (0\T{ VN (x)j K (y)f]e(0)}\a) , (1) 
where rj^, ?V and 7/0 are the interpolating fields associated with neutron, kaon and 9, 



respectively 



Hi- 



We next consider the expression (|1J) in terms of hadronic degrees of freedom and write 
the phenomenological side of the sum rule. Treating the kaon as a pseudoscalar particle, the 
interaction between the three hadrons is described by the following Lagrangian density: 

£ = igenK&JsKn for P = + 

£ = igenK®Kn for P = — (2) 

Writing the correlation function (0) in momentum space and inserting complete sets of 
hadronic states we obtain an expression which depends on the following matrix elements: 

-iV(p,p) = <n(p',s')\Q(p,s)K(q)>, 
(0\r) N \n(p',s')) = X N u s '(p') , 

(K(q)\j K \0) = Ajf, 
(e(p,a)\fj e \0) = X e u s (p) for P = + 

(e(p,s)\fj e \0) = -X e u s (p) l5 for P = - (3) 
Using the simple Feynman rules derived from (J2J) we can rewrite V(p,p') as 

V(p,p') = -gen K u s '(p'h 5 u s (p) P = + 

V(p,p') = -genKU s '(p')u s (p) P = - (4) 

The coupling constants and Ae can be determined from the QCD sum rules of the 
corresponding two-point functions. Xk is related to the kaon decay constant through 

X K = JH^. (5) 
m u + m s 



Combining the expressions above we arrive at 



—genK^e^N^K 
phen (p' 2 - m 2 N ) (q 2 - m? K+ ) (p 2 - m| ) 

x T E + continuum (6) 



with 



T E = a^^q^ - im N fa 
+ i(m N =f m e ) p"j5 

+ hb(p' 2 =F m e m N - qp') (7) 



We shall work with the cr^^q^pl structure because, as it was shown in [12[, this structure 
gives results which are less sensitive to the coupling scheme on the phenomenological side, 
i.e., to the choice of a pseudoscalar or pseudovector coupling between the kaon and the 
baryons. 

We now come back to and write the interpolating fields in terms of quark degrees of 
freedom as 

3k{v) = s(y)nMy): 

VN (x) = e abc (d T a (x)C ltl d b (x)h^ c (x) , 
^(0) = - e abc e def e cf9 s^(0)C 

x [J e (0) 75 C^(0)]K(0)C^(0)]. (8) 

The pentaquark current above (proposed in (|) contains a pseudoscalar and a scalar 
diquark. With these diquarks the two point function might receive a significant contribution 
from instantons. In [l3| we have studied a situation in which these instanton contributions 
affected the t wo-p oint function but gave a negligible contribution to the three-point function. 
Moreover, in jl4| we have observed that instantons give a negligible contribution to heavy 
baryon weak decays. Motivated by these results, in this first calculation we shall neglect 
instantons. 

Inserting the currents into (Q), the resulting expression involves the quark propagator 
in the presence of quark and gluon condensates, which is known from previous studies. 
Using it we arrive at a final complicated expression for the correlator, which is represented 
schematically by the sum of the diagrams of Fig. 1. 



Let us consider the phenomenological side (JdJ) and, following |15l |. rewrite it generically 

as: 

T(q 2 ,p 2 ,p' 2 ) = T pp + T pcl + T pc2 + T cc (9) 
where T pp (q 2 ,p 2 ,p' 2 ) stands for the pole-pole part and reads 

PP ip 2 — rn 2 e ){p' 2 — m 2 N )(q 2 — m 2 K ) 
The continuum-continuum term T cc can be obtained as usual, with the assumption of quark- 

hadron duality 

The pole-continuum transition terms are contained in r pcl and T pc2 . They can be explic- 
itly written double dispersion integral: 

bi(u,p 2 ) du 



^ pcl (m 2 N - p' 2 )(u - q 2 ^ 



, b 2 (s,p 2 )ds 
pc2 I ' (m 2 K -q 2 )(s-p' 2 ) ' { ' 



.2 

■ N * 



Since there is no theoretical tool to calculate the unknown functions bi(u,p 2 ) and &2(s,p 2 ) 
explicitly, one has to employ a parametrization for these terms. We will use two different 
parametrizations: one with a continuous function for the G and one where the pole term is 
singled out. 

We shall assume that the functions b\ and b 2 have the following form: 



~ poo 

bi(u,p 2 ) = b\{u) I duo 

Jrrik 



m| UJ — p 2 

b 2 (s, P 2 ) = b~( s ) r^-^H ( i2 ) 

Jm% UJ —p l 

with continuous functions bi >2 (w), starting from tuq. This is our parametrization A. 
The functions b\(u) and b 2 (s) describe the excitation spectra of the kaon and the nucleon, 
respectively. After Borel transform, the pole-continuum term contains one unknown constant 
factor which can be determined from the sum rules. 

In order to investigate the role played by the continuum, we shall now explicitly force 
the phenomenological side to contain only the pole part of the 0, both in the pole-pole term 
and in the pole-continuum terms. This can formally be done by choosing bi(u) = b 2 (uj) = 
8{uj — m@) in (H2J) and the functions then read: 



This is our parametrization B. In this case we have the in the ground state. Again, in 
the final expressions this gives additional constants which can be calculated. 



III. SUM RULES 

The sum rule may be written identifying the phenomenological and theoretical descrip- 
tions of the correlation function. As mentioned above, we shall work with the cr fiu 'j5q^p' u 
structure. In the case of the three-point function considered here, there are two independent 
momenta and we may perform either a single or a double Borel transform. We first consider 
the choice: 

(I) q 2 = p 2 = p' 2 (14) 

and perform a single Borel transform: p 2 = —P 2 and P 2 —>■ M 2 . In this case we take 
m 2 K ~ and single out the l/g 2 -terms. The second choice is: 

(II) q 2 ^0 p 2 =p' 2 . (15) 

Here we perform two Borel transforms: p 2 = —P 2 and P 2 —>■ M 2 and also q 2 = —Q 2 and 
Q 2 — > M' 2 . We have also considered the choice q 2 = p 2 = p' 2 = — P 2 , performing one 
single Borel transform (P 2 — > M 2 ). However, in the present calculation we were not able to 
find a stable sum rule. Introducing the notation G = —genK^e^N^K and using (I) and (II) 
we obtain the following sum rules: 
Method I: 

T PP (M 2 ) + T pc2 (M 2 ) = f ° dsp th (s) e- s ' M2 (16) 

Jo 

with 

T PP (M 2 ) = G - 2 5 (17) 

and for the pole-continuum part we obtain 

r pc2 (M 2 ) = Ae-<* /M2 param. A 

r pc2 (M 2 ) = Ae" m e/ M2 param. B (18) 

In both parametrizations the term Y pc i is exponentially suppressed and, as discussed in jl^ , 
has been neglected. A is an unknown constant and can be determined from the sum rules. 
Method II 

T pp (M 2 ,M' 2 )+T pc2 (M 2 ,M' 2 )= (19) 



H du f ° ds p th (s, u) e~ 3 ' M2 e- u ' M ' 2 (20) 
Jo Jo 



with 

-m|/M 2 _ -m%/M 2 



T pp = Ge-y™' 2 6 - 2 (21) 
and with 

T pc2 = A e -</ M ' 2 e"<* / m2 (22) 

for parametrization A and 

r pc2 = Ae-< /M ' 2 e- m e/ M2 (23) 

for parametrization B. Also in this case r pcl is exponentially suppressed. In the above 
expressions p t h is the double discontinuity computed directly from the theoretical (OPE) 
description of the correlation function (see [ill for details and also [16]) and so is the con- 
tinuum threshold of the nucleon defined as s = {m N + A N ) 2 . 



IV. RESULTS 



The hadronic masses are = 938 MeV, m^* = 1440 MeV, mx = 493 MeV and 
rriQ = 1540 MeV. For each of the sum rules above (Eqs. (|16|) and ([19))) we can take the 
derivative with respect to 1/M 2 and in this way obtain a second sum rule. In each case we 
have thus a system of two equations and two unknowns (G and A) which can then be easily 
solved. 

In the numerical analysis of the sum rules we use the following values for the condensates: 
(qq) = -(0.23 ± 0.02) 3 GeV 3 , (ss) = 0.8 (qq), < sg s a.Gs >= m 2 (ss) with m 2 = 0.8 GeV 2 
and (g 2 G 2 ) = 0.5 GeV 4 . The gluon condensate has a large error of about a factor 2, but its 
influence on the analysis is relatively small. The couplings constants Xn and Ae are taken 
from the corresponding two-point functions: 

X N = (2.4 ± 0.2) x 10" 2 GeV 3 (24) 

A e = (2.4 ±0.3) x lO^GeV 6 (25) 

The coupling A^ is obtained from (JHJ) with fx = 160 MeV, m s = 100 MeV and m u = 5 
MeV: 

X K = 0.37 GeV 2 . (26) 
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In Fig. 1, among these OPE diagrams there are two distinct subsets. In the first (from 
la to lg) there is no gluon line connecting the petals and therefore no color exchange. A 
diagram of this type we call color-disconnected. In the second subset of diagrams (lh, li 
and lj) we have color exchange. If there is no color exchange, the final state containing 
two color singlets was already present in the initial state, before the decay, as noticed in j^] . 
In this case the pentaquark had a component similar to a K — n molecule. In the second 
case the pentaquark was a genuine 5-quark state with a non-trivial color structure. We 
may call this type of diagram a color-connected (CC) one. In our analysis we write sum 
rules for both cases: all diagrams and only color-connected. The former case is standard 
in QCDSR calculations and therefore we shall omit details and present only the results. 
The latter case implies that the pentaquark is a genuine 5-quark state and the evaluation 
of genK will thus be based only on the CC diagrams. We shall work in the Borel window 
given by 1 GeV 2 < M 2 ,M' 2 < 1.5 GeV 2 . Since the strange mass is small, the dominating 
diagram is Fig. lb of dimension three with one quark condensate. In the range considered, 
the dimension 5 condensates are substantially suppressed compared to this term. 

We have found out that the contribution from the pole-continuum part is of a similar size 
as the pole part. For lower values of M 2 around 1 GeV 2 , the pole contribution dominates, 
however, for larger values of M 2 the importance of the pole-continuum contribution grows 
and eventually becomes larger than the pole part. This is an additional reason to restrict 
the analysis to small values for the Borel parameters. 

We have evaluated the sum rules for the coupling constant computed with all diagrams 
of Fig. 1 and we have found that they are very stable. We give the values of the coupling 
extracted at M 2 = 1.5 GeV 2 and M' 2 = 1 GeV 2 in Table I. In what follows we shall present 
our results for the coupling constant ge n K obtained with the color connected diagrams only. 
In Fig. 2 we show the coupling, given by the solution of the sum rule I A (I16|) . as a function of 
the Borel mass squared M 2 . Different lines show different values of the continuum threshold 
Ayv. As it can be seen, genK is remarkably stable with respect to variations both in M 2 and 
in An. In Fig. 3 we show the coupling obtained with the sum rule II A ()19j) . We find again 
fairly stable results which are very weakly dependent on the continuum threshold. In Fig. 
4 we show the results of the sum rule I B. In Fig. 5 we present the result of the sum rule II 
B. The meaning of the different lines is the same as in the previous figures. The results are 
similar to the cases before. 



S 



In Table I we present a summary of our results for ge n K giving emphasis to the difference 
between the results obtained with all diagrams and with only the color-connected ones. For 
the continuum thresholds we have employed A N = A K = 0.5 GeV. 



Table I:ge n K for various cases 



case 


fi'enx (CC) 


\ge nK \ (all diagrams) 


I A 


0.71 


2.59 


II A 


0.82 


3.59 


I B 


0.84 


3.24 


II B 


0.96 


4.48 



For our final value of genK we take an average of the sum rules I A - II B. It is interesting 
to observe that the influence of the continuum threshold is relatively small, especially when 
compared to the corresponding two-point functions. 

Considering the uncertainties in the continuum thresholds, in the coupling constants 
^k,n,0 and in the quark condensate we get an uncertainty of about 50%. Our final result 
then reads: 

\genK\ (all diagrams) = 3.48 ±1.8, 

|<7en*r|(CC) = 0.83 ±0.42. (27) 

Including all diagrams, the prediction for T@ is then 13 MeV (652 MeV) for a positive 
(negative) parity pentaquark. In the CC case we get a width of 0.75 MeV (37 MeV) for 
a positive (negative) parity pentaquark. The measured upper limit of the width is around 
5-10 MeV both in the Kn channel (considered here) and in the Kp channel. 

We see that it is very difficult to obtain the measured decay width for a negative parity 
pentaquark. 

V. SUMMARY AND CONCLUSIONS 

We have presented a QCD sum rule study of the decay of the 6 + pentaquark using a 
diquark-diquark-antiquark scheme with one scalar and one pseudoscalar diquark. Based on 
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the evaluation of the relevant three-point function, we have computed the coupling constant 
genK- In the operator product expansion we have included all diagrams up to dimension 5. 
In this particular type of sum rule a complication arises from the pole- continuum transitions 
which are not exponentially suppressed after Borel transformation and must be explicitly 
included. The analysis was made for two different pole-continuum parametrizations and in 
two different evaluation schemes. The results are consistent with each other. In addition, we 
have tested the idea8 preset in QQ by including o„. y diag^ with color exchange. 
Our final results are given in eq. (|27jl . 

We find that for a positive parity pentaquark a width much smaller than 10 MeV would 
indicate a pentaquark which contains no color-singlet meson-baryon contribution. For a 
negative parity pentaquark, even under the assumption that it is a genuine 5-quark state, 
we can not explain the observed narrow width of the 0. 
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FIG. 1: The main diagrams which contribute to the theoretical side of the sum rule in the relevant 
structure, a) - g) are the color-disconnected diagrams, whereas h) - j) are the color-connected 
diagrams. The cross indicates the insertion of the strange mass. 
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FIG. 2: l^enA - ! i n case I A with three different continuum threshold parameters. Solid line: 
A N = 0.5 GeV, dotted line:Ajv = 0.4 GeV, dash-dotted line: A^ = 0.6 GeV. 
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FIG. 3: \genK\ in case II A. Solid line: A N = 0.5 GeV. Dotted line: A N = 0.4 GeV. Dashed line: 
Atv = 0.6 GeV. M' 2 = 1 GeV 2 . 




FIG. 4: |<7enid i n case I B with three different continuum threshold parameters. Solid line: A^r 
0.5 GeV, dotted lineiAjv = 0.4 GeV, dash-dotted line: A N = 0.6 GeV. 
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FIG. 5: \genK\ in case II B. Solid line: A N = 0.5 GeV. Dotted line: A N = 0.4 GeV. Dashed line: 
Atv = 0.6 GeV. M' 2 = 1 GeV 2 . 
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